Charge order and phonon renormalizations: Possible implications for cobaltates. 
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Several experimental and theoretical studies in cobaltates suggest the proximity of the system to 
charge ordering (CO). We show, qualitatively, in the frame of a t — V model coupled to phonons 
that optical phonon modes at the K and M points of the Brillouin zone, which involves only O-ions 
displacement around a Co-ion, are good candidates to display anomalies due to the CO proximity. 
If by increasing of H2O content the system is pushed closer to CO, the mentioned phonon modes 
should show softening and broadening. 
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The discovery of superconductivity in Na x CoC>2 -yH-iO 
for x = 0.35 and y = 1.3 (Ref.fl) has motivated 
a large interest in the solid state physics community 
due, in part, to the similarities between this compound 
and high-T c cuprates. Cobaltates are 3d-electron sys- 
tems having a quasi- two-dimensional structure, C0O2, 
where the Co-atoms are in a triangular lattice. Re- 
cently, for describing superconductivity in these mate- 
rials, some authors proposed the interplay between elec- 
tronic correlations and the proximity to a charge den- 
sity wave (CDW) instability^^. These papers focus 
on a two-dimensional t — V model on a triangular lat- 
tice where the nearest neighbors Coulomb interaction 
V is the relevant parameter for bringing the system to 
the CDW instability which occurs for V > V c at the 
momentum q = (47r/3, 0). This CDW phase is called 
\/3 x \/3-CDW 2 ' 3 . Although some experimental reports 
in hydrated^ and unhydrated samples£&2 suggest a pos- 
sible interpretation in favor of the CO proximity, this 
scenario is still controversial^. Therefore, it is relevant 
to solve the CO controversy. In this paper we propose 
that if the system is close to the CO, signatures of it 
in the phonon subsystem can be expected. Hence, the- 
oretical and experimental studies on phonon anomalies 
may be important for testing the CO hypothesis. Those 
phonons which couple more directly with charge fluctua- 
tion are candidates to show anomalies and are discussed 
in this paper. 

Since the -\/3 x v^-CDW order takes place in the C0O2 
planes, we have performed for simplicity a 2D lattice dy- 
namics calculation^ in order to understand qualitatively 
which phonon modes are more sensitive to the proxim- 
ity of CO. For this purpose, we have studied the phonon 
eigenvectors. For q = we have obtained a double de- 
generated optical phonon, where only the O-ions move, 
associated with the Raman-active-in-plane E\ g mode-- 
of frequency u> ~ 480 cm -1 and also listed in table I 
of Ref. [HI . Increasing q along this branch the degener- 
acy is removed. At the iGpoint of the Brillouin zone 
[q = (47r/3,0)] we have obtained two phonon modes, 
A and B, with complex eigenvectors which lead to the 
displacements of the O-ions shown in FigQJ,. At the Ap- 
point [q = (-7T, 7r/v3)] there are two modes, C and D, with 
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FIG. 1: a) O-ions pattern of displacements of the two optical 
phonon modes at _ft"-point discussed in the text, b) The same 
for the M-point. Big (small full) circles represent the Co (O) 
ions. The crosses are the equilibrium positions of the O's. The 
arrows indicate the displacements of the O-ions. In the modes 
A and C the average distance to the Co-ion of the six O-ions 
changes with time leading to the possibility for the coupling 
of these modes to the density at the Co site. For B and D 
modes the average distance remains almost constant. The 
displacements of the O-ions in the mode A have a breathing- 
like pattern of displacements. 



real eigenvectors with O-ions displacements presented in 
FigOJ). Interestingly, in the modes A and C the aver- 
age distance to the Oo-ion of the six O-ions changes with 
time leading to the possibility for the coupling of these 
modes to the density at the Co site. In contrast, for B 
and D modes the average distance remains almost con- 
stant. As it is sketched in Figfjji the displacements of 
the O-ions in the mode A have a breathing-like pattern 
of displacements around Co-ions and therefore, an ade- 
quated symmetry to couple with the charge fluctuations 
at the Co-site when the system approaches the y/3 x \/3- 
CDW phase. In addition, this mode has the same q = K 
vector than the CO. 

For a qualitative study of those phonons which couple 
more directly with density fluctuations at Co site (A and 
C) we proposed the following HamiltoniarA 
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H = -t (^ itT Cj a + h.c.) + V ^ rnnj + ^^0(0^0^ + -) + g^(a^i + ai)rii (1) 
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FIG. 2: a) Phonon self-energy diagram. In II, the bare e-ph 
vertex, g (solid circle), is renormalized by electronic correla- 
tions as showed in b). Dotted and solid lines represent the 
propagators for phonons and fermions respectively. Double 
dashed line is the matrix-boson propagator- whose compo- 
nents Drb. and D\r are associated with the charge fluctua- 
tions and the enforcement of no-double occupancy constraint 
respectively. 

where t and V are the hopping and Coulomb repulsion 
between the nearest-neighbors sites i and j on the tri- 
angular lattice. cf iCT and Ci a are the fermionic creation 
and destruction operators for holes, respectively, under 



where the renormalized e-ph vertex 7(k, k — q; ioj n ) is 
7(k, k - q; iw n ) = -[D RR (q-uj n ) ^~ 

+^A_R(q; iu n ]\ (3) 

In ([2| Ao = 2g 2 N(0)/ujo is the bare dimensionless e- 
ph coupling and iV(0) is the bare electronic density of 
states (DOS). Ek = — tx(cos k x + 2 cos ^-ky cos ^k x ) is 
the renormalized electronic dispersion where x is the elec- 
tron doping. In the following t is considered to be 1 and 
energies are in units of t. We will present results for the 
commensurate doping x = 1/3 close to the value where 
the highest T c is obtained in cobaltates. For x = 1/3 
we obtain, in agreement with ARPES experiment o 17 ' 18 , 
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the constraint that double occupancy is excluded, n, is 
the fermionic density at the Co site. af i and are the 
phonon creation and destruction operators, respectively. 
g is the e-ph coupling for a given phonon of frequency 

UJQ. 

Theoretical calculations in Ref.[14( show that a one- 
band t — J model plus phonons can be derived, from the 
three-band Hubbard model plus phonons, for studying 
cases where the O-ions displacements lead to a modula- 
tion of the Zang-Rice singlet energy. This view was also 
used in Ref. 15] for studying the anomalous softening and 
broadening for the half-breathing mode in cuprates. 

The phonon self-energy corrections II(q, iw„) are com- 
puted after evaluating the diagram plotted in FigfS^,. The 
dark vertex in Figf2Ji is the e-ph vertex renormalized by 
correlations (Figf2jD). For the evaluation of the vertex 
we use a recently developed 1/iV-expansion for Hubbard 
operators^, where in order to get a finite theory in the 
N — > 00 limit, we rescaled t to t/N, V to V/N, and g to 
g/V~N and the spin projections a are extended from 2 to 
N. See Refs.JlU for details of the method for the present 
case. Note that phonon renormalizations are 0(1) which 
are of the same order than the free phonon propagator 
(Eq. (18) of Ref. 0]). 

Using Feynman rules and definitions of Ref. Q we ob- 
tain for the phonon self-energy n(q, iu) n ) 
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a large Fermi surface (FS) enclosing the T point. In a re- 
cent ARPES experiment 9 ., this FS topology was relevant 
for the interpretation in terms of the proximity of the 
system to the -s/3 x -s/3-CDW order. From we obtain 
for the dressed electronic DOS N*(0) - 0.631/t while the 
bare DOS is N(0) = x/2N*(0) - Q.ll/i. For the present 
parameters the instability to the v3 x V3-CDW order 
takes place for V > V c = 1.1-. 

The phonon frequency shifting and phonon linewidth 
are Aw = i?en(q, ujq + irf) and T = Imn(q, uj + irj) re- 
spectively. In Fig [3] we show Aw/Ao and T/Ao as a func- 
tion of V for the phonons at q = K and q — M . The fre- 
quency is ujq = t/2.5; being t ~ 150 meV , ujq ~ 480 cmT x 
which is close to the frequency of the two optical A 
and C phonon modes. With increasing V, approaching 
the V3 x y3-CDW phase, the two phonons become soft 
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(FigUi) and broad (FigH>). Since the x %/3-CDW 
instability occurs at q = K, softening and broadening of 
the q = K phonon are larger than for the q = M mode. 

The origin of the softening and broadening near the 
charge instability is understood as follows. The role 
played by charge densities is clearly evident by looking 
at eq.(3): As stated in Refs.JlQl Drr is proportional 
to the charge-charge correlation function % c (q, = 
— N(^) 2 Djui(q, u>) showing that charge fluctuations have 
a direct influence on the renormalized e-ph interaction. 
In FigJH are presented results for the density response, 
Jmx c (q = K,lo), for different V. The momentum was 
fixed at K where the V3 X V3-CDW phase takes place. 
For V = 0, it is clearly seen a collective peak at u> ~ A.2t 
at the top of the particle-hole continuum. With increas- 
ing V the collective peak softens and for V = 1 is still 
well defined at high energy of the order ~ 3t. But it 
has transferred spectral weight to a broad structure at 
low energies with a maximum at a scale ~ t/2. When 
the instability takes place, at V = V c = 1.1, this broad 
structure reaches u> = and the static charge suscepti- 
bility diverges. As result of this charge dynamics, the 
e-ph vertex is renormalized leading to the effects shown 
in FigO Since the instability is of static character, the 
effect will be stronger for low energy phonons than for 
large energy phonons. Contrary to the case for quarter 
filling on the square lattice^, the charge instability at 
x ~ 1/3 can not be seen strictly as the softening of a 
collective charge mode, but a redistribution of spectral 
weight takes place. The reason for this behavior can be 
found in the form of the particle-hole continuum which 
starts at u> = (inset in FigHJ and then, the collective 
peak never merges from the bottom of the continuum. 
Finally, it is important to notice that the static CO was 
not experimentally confirmed (see Reffioh. however the 
effects discussed in present paper need only a dynami- 
cal CO proximity in which the electronic system remains 
homogeneous. 

To our knowledge the exact value of the e-ph coupling 
is not known for cobaltates, however recent experiments 
suggest that this is nonegligeable^. In Refs.[l|3], for es- 
timations of T c , we have proposed for the total e-ph cou- 
pling A = 0.4 which is close to the order of other recent 
estimate d 21 ' 22 . In the context of superconductivity, this 
value must be related with A = A„ where v runs over 
the number of phonon branches and \ v is the coupling for 
each branch. Considering for instance Ao = 0.04, which 
is one tenth of A, we have found, for V ~ 0.9, a phonon 
frequency shifting Auj ~ 50 cm" 1 (15 cm" 1 ) for the K 
(M) phonon with respect to the situation for V << V c . 
A similar estimate can be done for T showing a broaden- 
ing of the order of T ~ 25 cm" 1 (6 cm" 1 ) for the K (M) 
phonon mode with respect to weakly or unhydrated sam- 
ples. In spite of the crudeness of the estimate, it shows 
that the effect is large enough as to be observed. 

It is important to make one remark. For q = 0, in the 
Ei g mode O-ions move parallel to the planes compressing 
the Co-atoms. However, from eq.(2) we expect phonon 
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FIG. 3: a) Phonon shifting Alu/Xq as a function of V for the 
phonon modes at q = K and q = M. The phonon frequency 
is ujo = t/2.5. b) Phonon broadening T/Xo as a function of 
V. With increasing V the system approaches the CO and the 
two phonons soften and broaden. 




FIG. 4: Imaginary part of the density response X c ( ( l = 
for x — 1/3, at K — (4tt/3, 0), for different V, approaching 
V c showing the redistribution of the spectral weight. For V 
near V c , low energy charge fluctuations renormalize the e-ph 
interaction and cause the phonon renormalizations shown in 
FigO Inset: particle-hole continuum. 



anomalies to be weaker for q ^ modes than for large 
q ones. Therefore, if q ~ phonons probe the charge 
fluctuations due to the proximity to CO, it is through 
inter-band transitions while, for large q phonons, intra- 
band transitions are directly involved and we expect the 
effect to be larger. 

Let us discuss how the presented results may be veri- 
fied experimentally. When studing superconductivity^ 
we obtain an increment of T c when V increases to V c . 
Based on the experimental results showing that T c in- 
creases with the water content y2£ we suggest a positive 
correlation between V and y. On the other hand, and in 
spite of a microscopic relation between V and y is lack- 
ing, experiments 2 ^, and theory 2 ^ show that hydratation 
causes the electronic structure to become more two di- 
mensional and then, one may argue that the Coulomb 
interaction between nearest-neighbors is less screened. 
Then, if hydration pushes the system closer to CO, we 
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expect larger phonon anomalies for samples with large y. 
Since the predictions in present paper are for large q op- 
tical phonons, they are in general accessible to neutron 
scattering experiments and, more recently, to inelastic 
x-rays scattering 22 . Since hydratation expands the lat- 
tice, mainly in the c-directior* 2 ^, one could expect that 
phonon dispersions also change due to this effect, making 
hard the detection of the predictions of the present pa- 
per. However, as the phonons modes discussed here are 
mainly of two dimensional character and perpendicular 
to the c-axis, we expect our result to be dominant. 

In summary, we have studied the influence of charge 
fluctuations, near the CO instability, on optical phonons 



modes at K and M points of Brillouin zone. We have 
discussed two in-plane optical phonons where only O- 
ions move compressing the Co- atoms. These phonons 
were predicted to be soft and broad when the systems is 
pushed near the charge instability. Softening and broad- 
ening were found to be larger for K than for M phonon 
because the former has the correct q and symmetry to 
couple with the y/3 x V3-CDW phase. 

The authors thank to M. Cardona, P. Lemmens, H. 
Parent, and M. Stachiotti for valuable discussions. We 
specially thank to A. Muramatsu for illuminating discus- 
sions. 
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